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Abstract 

We provide a detailed treatment of Ruijsenaars-Toda (RT) hierar- 
chy with special emphasis on its the theta function representation of all 
algebro-geometric solutions. The basic tools involve hyperelliptic curve 
1C P associated with the Burchnall-Chaundy polynomial, Dubrovin-type 
equations for auxiliary divisors and associated trace formulas. With 
the help of a foundamental meromorphic function </>, Bakcr-Akhiezcr 
vector $ on /C p , the complex- valued algebro-geometric solutions of RT 
hierarchy are derived. 

1 Introduction 

Nonlinear integrable lattice systems have been studied extensively in 
relation with various aspects and they usually possess rich mathematical 
structure such as Lax pairs, Hamilton structure, conservation law, etc. The 
Toda lattice is one of the most important integrable systems [10] . It is 
well-known soliton equations such as the KdV, modified KdV, and nonlin- 
ear Schrodinger equations are closely related to or derived from the Toda 
equation by suitable limiting procedures [TOj, [12]. Various kinds of Toda 
lattice have been discussed since it was proposed [8] EJ HUJ QT] H5] . Among 
them, a remarkable discovery was made by Ruijsenaars in the area of inte- 
grable lattice systems [T]. He found a relativistic integrable generalization 
of non-relativistic Toda lattice through solving a relativistic version of the 
Calogero-Moser system. The Lax representation, inverse scattering problem 
of the Ruijsenaars-Toda lattice and its connection with soliton dynamics 
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were investigated. A general approach to construct relativistic generaliza- 
tions of integrable lattice systems, applicable to the whole lattice KP hier- 
archy, was proposed by Gibbons and Kupershmidt [6]. After that, a series 
of techniques to construct relativistic lattice equations were developed by 
systematic procedure and Hirota's bilinear method [31 [T7] . 

The Ruijsenaars-Toda lattice, sometimes also called relativistic Toda 
(RT) lattice, takes the form [3] 

Pt = (l + hP)(a-ar), 

at = a (/3 + — /3 + fra + — ha~) 

in Flaschka variables or 

e x k+1 -x k 

x ktU = (i + hx k>t ) (i + hx k+ljt ) h2pXk+1 - Xk 
- (i + hx k ^ t ){\ + nx^)^-^—-— 

in Newtonian form, where the small time step h = c" 1 and c is light speed. 
In the non-relativistic limit c — > oo, the RT equation (jl.ip reduced to the 
well-known Toda lattice equation [7], 

fit = a — a~ , at = a((3 + — (3) (1-3) 

in Flaschka variables, or 

x kn = e x k+i-x k _ e x fe -x fe -i ( L4 ) 

in Newtonian form. Eq (jl.ip is the Poincare-invariant generalizations of the 
Galilei-invariant Toda systems (jl.3p . 

Mathematical frame work such as Lax representation, Backlund transfor- 
mation, Hamiltonian structure of RT lattice Eq (jl.ip . etc, were investigated 
by some authors |13j-|16j. Cosentino obtained a soliton solution by using the 
1ST method. Hietarinta and Junkichi Satsuma transformed the RT eq (jl.ip 
into trilinear form through a suitable dependent variable transform. Later 
Yasuhiro Ohta, etc. decomposed the RT lattice eq (jl.ip into three Toda 
systems, the Toda lattice itself, Backlund transformation of Toda lattice, 
and discrete time Toda lattice and explicitly derived the solutions in terms 
of the Casorati determinant [17] . The solution they obtained converges to 
that of TL eq (jl.3p in the limit of c — > oo. 

Algebro-geometric solutions (finite-gap solutions or quasi-period solu- 
tions), as an important character of integrable system, is a kind of explicit 
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solutions closely related to the inverse spectral theory [271129]. Around 1975, 
several independent groups in UUSR and USA, namely, Novikov, Dubrovin 
and Krichever in Moscow, Matveev and Its in Leningrad, Lax, McKean, van 
Moerbeke and M. Kac in New York, and Marchenko, Kotlyarov and Kozel in 
Kharkov, developed the so-called finite finite-gap theory of nonlinear KdV 
equation based on the works of Drach, Burchnall and Chaunchy, and Baker 
[25| l28j [30] . The algebro-geometric method they established allowed us to 
find an important class of exact solutions to the soliton equations. As a 
degenerated case of this solutions, the multisoliton solutions and elliptic 
functions may be obtained [271 EH] • Its and Matveev first derived explicit 
expression of the quasi-period solution of KdV equation in 1975 |28j . which 
is closely related to the finite-gap spectrum of the associated differential 
operator. Further exciting results appeared later, including the finite-gap 
solutions of Toda lattice, the Kadomtsev-Petviashvili equation and others 
[3 [30j [35], which could be found in the wonderful work of Belokolos, et al 
[27] . In recent years, a systematic approach based on the nonlinearization 
technique of Lax pairs or the restricted flow technique to derive the algebro- 
geometric solutions of (1 + 1)— and (2 + 1)— dimensional soliton equations 
has been obtained [31]- [32]. An alternate systematic approach proposed by 
Gesztesy and Holden can be used to construct algebro-geometric solutions 
has been extended to the whole (1+1) dimensional continuous and discrete 
hierarchy models [36j-[38j. |40l 141] . 

In this paper, we mainly discussed the algebro-geometric quasi-period 
solutions of RT hierarchy for fixed constant h = c _1 ^ 0. In following sec- 
tion 2, the RT equation (jl.ip is extended to a whole RT hierarchy through 
the polynomial recursive relations. In section 3, we give a detailed study 
of algebro-geometric solutions for the stationary RT hierarchy. Firstly we 
derived the hyperelliptic curve K, p in connection with the stationary RT 
hierarchy. Then a fundamental meromorphic function (j) on fC p , the Baker- 
Akhiezer vector ^ and the common eigenfunction of zero-curvature pair 
U, V p , was introduced to study the trace formula and asymptotic properties 
of <f> and ip\ , respectively. With the help of Riemann theta function associ- 
ated with /C p , one finds the theta function representations for 4> and ipi by 
alluding to Riemann's vanishing theorem and the Riemann-Roch theorem. 
In section 4, we derive the complex-valued algebro-geometric solutions of 
RT hierarchy with a given inital value problem by using the results in sec- 
tions 3 and 4. Finally, in Appendix A and B we give Lagrange interpolation 
representation and asymptotic spectral parameter expansions that will be 
used in this paper. 
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2 The Ruijsenaars-Toda Hierarchy 

In this section, we derive the Ruijsenaars-Toda hierarchy by using a polyno- 
mial recursion formalism. Throughout this section let us make the following 
assumption. 

Hypothesis 2.1. In stationary case we assume that u : M — > C satisfies 

a,/3£C z , a(n)^0, n£Z (2.1) 

In the time-dependent case we suppose u : M 2 — > C satisfies 

a(-,t), p(-,t) G C z , a(n, t) / 0, (n, t) G Z x K 

a(n, •), /3(n, •) G C 1 (IR). (2.2) 

In this paper, we denote by S the shift operators acting on ip = 
{i/}(ri)}n=-oo G C z according to (5 ± V)( n ) = VK n ± 1); or V' ± = S^tp for 
convenience. 

We start from the following 2x2 matrix isospectral problem |26] 



S+^ = U^ i>=[ ^ ), U=( ° 1 ), (2.3) 

where the functions a, (3 are potentials, and z is a constant spectral param- 
eter independent of variable n. 

Define sequences {fe}i€N an d {gi}eeN recursively by 

go = -1/2, f = 0, (2.4) 

ft+i + Pft + 9t + 9i = 0, ^N , (2.5) 

f^ft+i ~ ha+ fe+i ~ a fi + a+ ft ~ P9e + P9t 

-9e +1 + 9e+i = 0, I G No- (2.6) 

Explicitly, one obtains 

ft = 1, gi = 

f 2 = -h{a + + a) - P + 6 U 

92 = -h 2 a ++ a + -h 2 (a + ) 2 + h 2 a + a-ha + P + -ha + P -a + 

+ ha + 5 1 -5 2 /2, 

f 3 = h 2 a ++ a + + h 2 (a + ) 2 -h 2 a + a + ha + (3 + + 2ha + p + a + (2.7) 

—I— OOO / _i_ 

+ h a a + h a — H aa~ + 2haP + haP~ + a — (ha 
+ Tia + /3)<5i + 82, etc. 
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Here {5e}eei>i denote summation constants which naturally arise when solv- 
ing ()2.4p . Subsequently, it will be useful to work with the corresponding 
homogeneous coefficients fe,gi defined by the vanishing of all summation 
constants 5 k for k = 1, ,£, 

9o = -1/2, gt = ^1^=0^=1,-/, t G N, 
fo = 0, f t = fe\ Sj =o,j=i,-,e, ^N. 
By induction one infers that 

i 

gt = ^ j bt- k g k , i G N 

s =° (2.8) 

s=0 

introducing 

«5o = 1. 

Remark 2.2. The constants 5j (j G No), can be expressed in terms of the 
branch points Ei, (i = 0,- • • ,p + 1) of the associated spectral curve defined 
in (12T3UD . (Theorem 6) 

In order to obtain the RT hierarchy associated with the spectral problem 
(|2.3p . we first solve the stationary zero-curvature equation 



UV P - V+U = 0, y p = 0^)2x2 (2.9) 

with 



V11 v 12 



V v=\ „_ Z > (2- 10 ) 



\ ^21 ^22 / 
where each entry is a polynomial in z, 

P+l p+l 
Vh = ^ 5p+1 _ J z J +/ p+2 , y 12 = (2.11) 
i=o j=o 

F 21 = (fe-l)a+y+ = (hz-l)a+ ^pjfz p+1 ~^j , (2-12) 
^22 = -Y,9 P+ i-jZ 3 - (2.13) 

j=0 
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Equation (|2.9p can be rewritten as 

V 2l - (hz - l)aV 12 = 0, (2.14) 

V--V 11 -(z + f3)V 12 = 0, (2.15) 

(hz - l)aV£ + (z + j3)V 2l - (hz - l)aV 22 = 0, (2.16) 

(hz - l)aV{~ 2 + (z + p)V 22 - V 2l -{z + p)V 22 = 0. (2.17) 

Since det(U)^ for z £ C\{l/h} by $T[$, QEJ) yields tr(V+) = t^UVpU" 1 ) - 
tr(V^) and hence 

^11+^22 = ^11+^22, 

implying V\\ + V 22 is a lattice constant. If Vn + V 22 = c ^ 0, one can add a 
polynomial times the identity to V p , which does not affect the zero-curvature 
equation. Therefore we can choose 

Vn + V 22 = (2.18) 

without loss of generality. This fact leads to the following result. 

Lemma 2.3. Suppose U and V p satisfy the stationary zero-curvature equa- 
tion P^j) . Then (pD^ - fpU?]) change into 

V{i + V n + (z + P)V 12 = 0, (2.19) 
(hz - l)aV{~ 2 - (z + 0)V£ - (hz - l)a + Vf 2 + (z + p)Vu = 0. 

(2.20) 

In particular, the coefficients {fe}e=o : ... jP and {ge}i=o,...,p defined in (|2.1ip - 
(|2.13p satisfy the recursive relations (|2.4p - (|2.6p . 

Proof. Eq. (^T9]) and eq. (j!T20|) arise from (pl^) - (pT7|) by substituting V 22 
for -Vn. Insertion of (j2TTT]) - (j2TT3Tl into ([2TT9]) and (l230]l then yields the 
relations ([13])- $2M). □ 
Inserting (|2.1ip - (l2.13p into (|2.9p and using the results of Lemma l2.3l yield 
the following theorem. 

Theorem 2.4. Suppose that U and V p satisfy the stationary zero-curvature 
equation (|2.9p . Then (j2.9|) reads 

o = uv p -v+u 



( 

(Tiz - 1) (a/p+2 -a/p+i + 
V -af P+2 ) -/3(g~ +l - g p+ i) 



(2.21) 



G 



which is equivalent to 



fp+2 ~ /p+2 - 0> 
(3(9p +1 ~ 9 P +i) + af~ +l 



0. 



(2.22) 
(2.23) 



Thus, varying p G No, equations (|2.22|) and (|2.23p give rise to the sta- 
tionary Ruijsenaars-Toda (RT) hierarchy which we introduce as follows 



s-RT p (a, (3) 



fp+2 ~ fp+2 



0, peN 



(2.24) 

In the special case p = 0, one obtains the stationary version of the Ruijsenaars- 
Toda system (jl.ip 



/ ha(a~ - o+) + a(p~ - 0) \ 
y fo/3(a — q + ) — q + a + J 

In the case p = 1, one finds 

f a(h 2 a ++ a+ + h 2 (a+) 2 - h 2 a+a + ha+ (3+ + a H 
+haf3 — ha(3~ — h 2 aa~ — h 2 (a~) 2 + h 2 a~a~~ 
—2ha~/3~ — ha~ (3 — a~) 



—ha 2 — haa — af3 + ha + a ++ + h{a + ) 2 

+a + f3 + + f3(-h 2 a + a - h 2 a 2 + h 2 aa~ 
-haft - haj3~ - a + h 2 a ++ a + + h 2 (a + ) 2 
-h 2 a + a + ha + /3+ + ha + (3 + a + ) 

a(-ha + + UoT + P~) 

a — a + + ha/3 — ha + (3 



(2.25) 



0. 



In accordance with our notation introduced in (|2.8p . the corresponding ho- 
mogeneous stationary Ruijsenaars-Toda equations are defined by 

s-RT p (q, (3) = s-RT p (a, (3) \s a =i,8 e =o,e=i,..., P , P G N . 

Next we turn to the time-dependent Ruijsenaars-Toda hierarchy. For 
that purpose the potentials a and f3 are now considered as functions of 
both the lattice point and time. For each equation in the hierarchy, that 
is, for each p € No, we introduce a deformation (time) parameter t p € R in 
a,/3, replacing a(n),(3(n) by a(n,tp),/3(n,tp). The quantities {fe}ie^o an d 
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{9i}ieNo are s * m defined by f)2.4[) - (|2.6j) . The time-dependent Ruijsenaars- 
Toda hierarchy then obtained by imposing the zero-curvature equations 



Ut p (t p ) + U(t p )V p (t p ) - V p + (t p )U(t p 
Relation (|2.26p implies 





0, t p € 



(2.26) 



+ 



(hz-l){af p+2 
-afp+2) 



' a f P +i+ a+ fp+i 
-P(9p + i~9p+i) 



0. (2.27) 



RT p (a,/3) 



I 

\ [hz - l)a tp p tp 

Varying p € No, the collection of evolution equations 

a tp ~ a{f p+2 + f~ +2 ) \ 

v Ptp ~ P(9p+i ~ 9p+i) ~ «/ P +i + « + / P +i / 

(n,t p ) £ZxR, p € No, (2.1 

then defines the time-dependent Ruijsenaars-Toda hierarchy. Explicitly, 

at Q — ha(a~ — a + ) — a(/3~ — j3) 
fit — hf3(a — a + ) + a — a + 



RT («,/3) 



0. 



RTi(a,/3) 



+ <5i 



a tl — a(h 2 a ++ a + + h 2 (a + ) 2 - h 2 a + a 
+ha + f3 + + a + + haft — haj3~ — h 2 aa~ 
-h 2 (aT) 2 + h 2 a~a - 2ha~j3~ 
—ha~P — a") 

fiti + ha 2 + haa~ + aj3~ — ha + a ++ — h(a + ) 2 
—a + (3 + — fi(—h 2 a + a — h 2 a 2 + h 2 aa~ 
-hafi - hafi- - a + h 2 a ++ a+ + h 2 {a + ) 2 
-h 2 a + a + ha + j3 + + ha + f3 + a+) 

-a(— ha + — (3 + ha~ + /3~ 



\ 



-(a — a + + hafi — ha + (3) 



0, etc., 



represent the first two equations of the time-dependent Ruijsenaars-Toda hi- 
erarchy. The system of equations, RTo(a, /3) = is of course the Ruijsenaars- 
Toda system ([l.lj) . 

Next, taking into account (|2.18p . one infers that the expression R 2p+2 , 
defined as 

R 2p+2 (z) = -V^{z,n) - Vi 2 (z,n)V 2 i(z,n), 

= -V^(z,n)-(hz-l)a + V 12 (z,n)V^(z,n) (2.29) 
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is a lattice constant, that is, i?2p+2 — -^2p+2 = 0, si nce taking determinants 
in the stationary zero-curvature equation (|2.9p immediately yields 

(Hz - l)a {{-V{[) 2 - V£Vji + V& + VvMn) = 0. 

Hence, i?2p+2(^) only depends on z, and one may write i?2p+2 as 

1 2p+l 

#2p+ 2 (*) = — J] ( Z " {ErnffiS 6 C, J> G N . (2.30) 

m=0 

Relations (|2.29p and (|2.30p allows one to introduce a hyperelliptic curve JC P 
of (arithmetic) genus p (possibly with a singular affine part), where 

2p+2 

Kp : T p {z, y) = y 2 + AR 2p+2 = y 2 - Y[(z - E m ) = 0, p G N . (2.31) 

m=0 

Equations (|2.19|) . (|2.2(jp and (|2.29p permit one to derive nonlinear dif- 
ference equations for Vn,Vi 2 separately. One obtains 

_V* - (hz - l)a+(V+ + Vu)(Vu + V n ) 

= {z + P)(z + P+)R 2p+2 , (2.32) 

(ft? - l)(a^2 - a + V+) + {z + f3) 2 V l2 \ 2 

-4(fe -l)(z + f3) 2 a + V 12 V 1 + 2 = 4(z + P) 2 R 2p+2 (2.33) 

Equations analogous to (I2.32P and (12.33H can be used to derive nonlinear 
recursion relations for homogenous coefficients fi,g£- In addition, as proven 
in Theorem B.l, (I2.32P leads to an explicit determination of the summa- 
tion constants 8i, S 2 , ■ ■ ■ , S p in (|2.24p in terms of the zeros Eq, . . . , E 2p +i of 
associated polynomial R2 P +2 in (|2.30p . In fact, one can prove (cf. (|4.92|) ) 

5 i = c e (E), £ = 0,...,p, (2.34) 

where 

2p+l 

2 



1 ^p+i 

co(K) = l, Cl (E) = --Y,En 



m=0 



Ck(E) = Yl 



(2jo)!...(2j 2 p +1 )!^ ---^ + + 1 1 
2 2fc (joD 2 • ■ ■ (i2p + i0 2 (2j o -!)••• (2j2p+i - 1) 



j0vJ2p+l=0, 
jo + "-+j2p+l=fc 

k G N. 

are symmetric functions of E = (Eq, . . . , E 2p +{). 
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Remark 2.5. If a, f3 satisfy one of the stationary Ruijsenaars-Toda equa- 
tions in (|2.24p for a particular value of p, s-AL p (a,/3) = 0, then they satisfy 
infinitely many such equations of order higher than p for certain choices of 
summation constants 5g. This is seen as follows. Assume f p +2 — f p+2 = 
for some p G N and some set of integration constants {5e}e = i t ,,, tP C C, one 
infers 

fp+2 = Ap+2 

for some constant A p+ 2 6 C Subtracting the constant X p +2 (i.e. writing 
fp+2 = J2^=o^p+2-kfk, for some set of constants {Se}e=i,..., p +2 and absorb- 
ing Xp+2 into §p+2), we may without loss of generality assume that f p+ 2 = 0, 
and hence the recursion (|2.6p implies 

P{g p+ i - 9 P +i) + otfp +1 - a + / p + +i = 
in (I2.24p is equivalent to 

9p+2 ~ 9 p+2 - ha+ f p +2 ~ Ka fp+2 = °- 

Hence, 

9p+2 = A* +2 

for some constant A* +2 G C. This indicates 

^9 p +2 ~ 9 P +2) + a/p+2 - « + /p+2 = °- 

Then 

fp+3 = Ap+3 

for some constant A p +3 6 C which arises from (|2.5p . Similarly, subtracting 
the constant A p+ 3 (i.e. writing f p+ ^ = YT s =o&p+z-kfki f° r some set of con- 
stants {^}^=i,...,p+3 and absorbing A p+ 3 into 6 p+ s), we may without loss of 
generality assume that / p+ 3 = 0, . . . Iterating this procedure yields 

s-RT q (a,/3) = 

for all q > p + 1 (corresponding to some p-dependent choice of integration 
constants {8t]e.=i,...,p)- 
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3 Stationary Algebro-geometric Solutions 



This section is devoted to a detailed study of the stationary Ruijsenaars- 
Toda hierarchy and its algebro-geometric solutions. Our basic tools are 
derived from combining the polynomial recursion formalism introduced in 
Section [2] and a fundamental meromorphic function on a hyperelliptic 
curve fC p . We will obtain explicit Riemann theta function representations 
for the meromorphic function <j), the Baker- Akhiezer function ipi, and the 
algebro-geometric solutions a, (3. 

Unless explicitly stated otherwise, we suppose in this section that 



and assume ([IlJ-flSID, (t^9l)- (f2TTD . ([2TT9]) - ([2T24l . (p^9l) - (p3TT) . keeping p 6 
N fixed. 

Throughout this section we assume K p defined in (12.311) to be nonsingu- 
lar, that is, we suppose that 

E m + E m , for m ^ m', E m e C\{h}, m = 0, . . . , 2p + 1. (3.2) 

We compactify IC p by adding two points ^00+ and Poo--, P00+ 7^ -P00-5 at 
infinity, still denoting its projective closure by K, p . Finite points P on K, p 
are denoted by P = (z, y) where y(P) denotes the meromorphic function on 
IC P satisfying F p {z,y) = 0. The complex structure on IC P is then defined in 
a standard manner and 1C P has topological genus p. Moreover, we use the 
involution 



*:K P ^K P , P=(z,y)^P* = (z,-y), P^ ^ P^ ± = P^. (3.3) 



We also emphasize that by fixing the curve IC P (i.e., by fixing Eq, . . . , £2p+i), 
the summation constants {^}£=o,...,p in the stationary RT p equations are 
uniquely determined as is clear from (|2.34[) . which establish the summation 
constants 5e as symmetric functions of Eq, . . . , E2 P +i- 

For notational simplicity we will usually tacitly assume that p € N.(The 
trivial case p = is explicitly treated in Example 3.6) 

In the following, the zeros of the polynomial V\2{-,n) (cf. (12. lip ) will 
play a special role. We denote them by {A t j( n )}j=i,.--,p an< ^ hence write 



(3.1) 



p 



v 12 (z)=n(z- H ). 



(3.4) 
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Similarly we write 

1 P 

V 2 i(z) = ha + (z--)l[(z-fip, / x+(n)= A t,(n + l), (3.5) 
i=i 

j = l,...,p, n G Z, 

and we recall that (cf. (f^29]) ) 

P 2p+2 + = -{hz - l)aV 12 V+. (3.6) 

The next step is crucial; it permits us to "lift" the zeros fij from the complex 
plane C to the curve /C p . From (|3.6p one infers that 

R2 P +2{z) + Vn(z) 2 = 0, z G {Mj,M/,lA}j=i,...,p- 

Now we introduce {/Uj}j = i i .... p C /C p , {^}j=i,..., p C /C p and P^ G /C p by 

/ij(n) = (pj(n),-2Vu(pj(n),ri)), j = 1, . . . ,p n G Z, (3.7) 
A J + W = K + W,2y 11 (/i+(n),n)), i = l,--.,p nGZ, (3.8) 



and 
where 



P ft = (l/M^n(l/£,n)), (3.9) 



2p+l 

(2y 11 (l// l ,n)) 2 =y(l/^) 2 = JJ(ft-E m ) 

m=0 

is independent on n. 

Next we briefly define some notations in connection with divisors on /C p 
[371 [38]. A map, V : K p -> Z, is called a divisor on /C p if V(P) ^ for only 
finitely many P G /Cp. The set of divisors on /C p is denoted by Div(/C p ). We 
shall employ the following (additive) notation for divisors, 

*DqoQ = + T>q, Vq = Vq 1 + . . . + V Qm , 

Q = {Qi, • • • , Qm} e Sym m /C p , Q G £ p , m G N, 



where for any Q G /C p , 

P Q :/C P ^N , P^V q (P) 



1 for P = Q 

for Pg/C p \{Q}, 



and Sym n /C p denotes the nth symmetric product of /C p . In particular, one 
can identify Sym m /C p with the set of nonnegative divisors < T> G Div(/C p ) 
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of degree m. Moreover, for a nonzero, meromorphic function / on K p , the 
divisor of / is denoted by (/). Two divisors V, £ € Div(/C p ) are called equiva- 
lent, denoted by V ~ £ , if and only if T> — £ = (/) for some / € A4(/C p )\{0}. 
The divisor class [V] of V is then given by [V] = {£ € Div(/C p )j£> ~ £}. We 
recall that 

deg((/)) = 0, f e M(K p )\{0}, 

where the degree deg(P) of V is given by deg(P) = J2petc p ^(-f )■ 

Next we introduce the fundamental meromorphic function on K p by 

= Vl2 (z,n) (3 - 10) 

= ihz -f' V r iz \ n \ P = (z,y)eK p ,neZ, (3.11) 
y/2 + Vh(z,n) 

with divisor (<^>(-,n)) of cf>(-,n) given by 

(0(-,n)) = P P ^ +(n) - ^.AH, (3.12) 
using (|3.4p and (|3.5p . Here we abbreviated 

A = A + = {At, • • • ,Ap }• 

Given <j)(-,n), the meromorphic stationary Baker-Akhiezer vector \I/(.,n, no) 
on IC P is then defined by 

,wd Z' ^i(^ n ' n o) \ fo ,a\ 

*(P,n,n ) = , (3.13) 



V>i(-P,n,n ) 



Y[ n n 7l n J{P,n% n'>n , 

1, n' = n , (3.14) 
lUXn^P^T 1 , n'<n , 

n:ri rto+1 (^#+. + /3(n0), n'>n , 

^ 2 (P,n,n ) =^(P,no) x <( 1, n' = n , 

EGSUh-i (^ffe 11 + ^ + fl"')) ~\ < no- 



(3.15) 

Basic properties of <^> and ^ are summarized in the following result. 



Lemma 3.1. Suppose a, j3 satisfy (|3.ip and the pth stationary Ruijsenaars- 
Toda system (1X241 . Moreover, assume (12T3UD (f23Tp and $FZ$) and let 
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P = (z,y) E }C P \{P 00 ±, Pfr}, (n,no) € 1? . Then cf> satisfies the Riccati-type 
equation 

as well as 



<f>(P)<fT(P) -{z + P)<t>~(P) - (hz - l)a = 0, 

(Hz - l)a+V&(z) 



ct>{p) + <j>{p*) 
0(P)-0(P*) 



^12(2) 



tt 2 (z) 



The vector $ satisfies 

ip2(P,n,n ) = ifji(P,n, n )(j)(P, n), 
U(z)W-(P) = *(P), 
y p (z)^(P) = (y/2)^(P), 

^i(P,n,n )Vi(P*,n,n ) = (ft? - l) n - n °r(a+ n, n ) 

Vi2(z,n) 
^12(2, n )' 

7pi(P,n,no)^2(P*,n,n ) + ipi(P* ,n,n )ip 2 (P,n,n ) 

V n (z,n) 



-2(hz-l) n - n °T(a + ,n,n ) 

) - tpi(P*,n } n )ip2(P,n, 
-(hz-l) n - no T(a+,n,n ) 



Vi2(z,n Q ) 

ipi(P,n,n )ip2(P*,n,no) - ipi(P*, n, n )^ 2 (P, n, n ) 



^12(2,^0) 



where we used the abbreviation 



r(/,n,n ) = < 



Un>=n f( n ')> n>n 0l 
1, n = n , 



(3.16) 

(3.17) 

(3.18) 
(3.19) 

(3.20) 
(3.21) 
(3.22) 

(3.23) 

(3.24) 
(3.25) 



(3.26) 



iro=n/(n')-\ n<n . 
Proof. To prove f|3. 16j) one uses the definitions (|3.10j) of <fi and equations 
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(12J9|) ([gjZOD and ([239]) to obtain 

4>{P)4T{P) ~(z + P)<F{P) ~ (hz - l)a 



y/2 - y/2 - F lT y /2 - V£ _ _ 

V12 V{~ 2 H) Vn 

1 

^12^12 

- (Hz - l)Vi2Vn 



(y/2 - V u )(y/2 - V{{) - (z + (3)V 12 (y/2 - V{[) 



= 0. 

Equations (|3.17p - (|3.19p are clear from the definitions of <j) and y. Next we 
use induction to prove (|3.20p . 
(i) n = uq; one easily finds 

tfj 2 (P,n ,n ) = <p(P,no) = tpi(P, n , n o )0(P, n ). 

by definition of 

(ii)n > no; we assume (|3.20p holds for n = no, . . . , n — 1. Then ip\,ip2 satisfy 

ip2(P,n,n ) _(hz — l)a«/)~(P, n) -1 + z + /3 ^(P, n, n ) 
V>i(P,n,n ) <P~(P,n) tp±(P,n,n ) 

= ,_, p v ^ (P,n,n ) 

(P,n) 

=(fe- l)Q^~(P,n)- 1 + z + (3.27) 



that is, 

■0 2 (P,n,n o ) 



0" (P, n) - (z + P)<fT (P, n) -{hz-l)a = 0. (3.28) 



ipi(P,n,n ) 

Comparing (|3T2"gj) with (f57TU|) then yields (f3T2U|) for n > n , n € N. 
(iii)n < no; analogous proof with (ii). 
The definition of eft (cf. (|3.14p ) implies 

MP,ri,n ) =^(P,n,n )4 ) -(P,n) (3.29) 

and hence 

tl)\{P,n,7iQ) = ^(P,n,n ), (3.30) 
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which follows from (|3.20p and (|3.29p . The definition of ^2 (cf. ()3.15[) ) implies 

/ (hz — l)a \ 
ip 2 (P, n, n ) = I ^zrrp ^ + z + pjip 2 (P, n, n ) 

= (hz - l)aipi (P, n,n ) + (z + /3)-0^"(P, n, n ), (3.31) 

where we use (|3,20p again. Then equation (|3.2ip follows from (|3,30p and 
p.3ip . Property (|3.22|) is an immediate consequence of (|3.20p and the def- 
inition of (j). Finally, Equations (l333l) - (^25|) follow from (l3T7D - (l3T9l) ,the 
definition of Vi (cf. (IBTTi]) ) and (^201) . □ 

Combining the polynomial recursion approach in the section [2] with 
(|3.4p yields the following trace formula, which means fi, gi can be expressed 
by the symmetric functions of the zeros fij of V± 2 . For simplicity, we only 
show one of them. 

Lemma 3.2. Suppose that a, f3 satisfy the pth. stationary Ruijscnaars-Toda 
system (f!T2"4"]) . Then, 

p 

-h(a + a + ) - P + 8 1 = -J^^r ( 3 - 32 ) 

i=i 

Proof. Relation (|3.32p are proved by comparison of powers of z equating 
the corresponding expression (]3.4|) for V± 2 with that in (|2.1ip and with (|2.7p 
taken into account. □ 

Next we turn to asymptotic properties of (f> and ipi in a neighborhood 
of Poo± and P^ The asymptotic behavior of $> 2 is derived naturally from 
(I3.20p . This is a crucial step to construct the stationary algebro-geometric 
solutions of stationary Ruijsenaars-Toda hierarchy. 

Lemma 3.3. Suppose that a, /3 satisfy the p-th stationary Ruijsenaars- 
Toda system (|2.24p . Moreover, let P = (z,y) G /C p \{Poo ± , P h }, (n,n ) € 
ZxZ. Then, (ft has the asymptotic behavior 

_ rr' + OJ-taj+oK), c = 1/z , 

C^o + (/i 2 a+a + a- fe/3)C + 0(C ), P -> -P00+, 

(3.33) 

,=„ ^T^ + °(C 2 ), P->fli, C = z-l/h, (3.34) 
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Accordingly, the component ip\ of the Baker-Akhiezer vector have the 
asymptotic behavior 



wi{P,n,nn) = < , C = l/z, 

c-o lr(fto+) (i + O(0), p^p*c + . 



(3.35) 

P^P ft , C = 2- V^(3.36) 



The divisor (ipi) of is given by 

(ipl(;n,n )) = £>£ (n) - V^ {no) + (n - n )(Dp, - V Poo _). (3.37) 

Proof. The existence of the asymptotic expansions of 4> m terms of the 
local coordinate Q = 1/ z near Poo±, respectively, Q = z—l/h near P^ is clear 
from the explicit form of (j) in (|3.10j) and (|3. 1 1 [) . Insertion of the polynomials 
(|2.1ip and (12.12P into (13.10P and (|3.1ip then yields the explicit expansions 
coefficients in (I3.33|) and (j3.34p . Alternatively, and more efficiently, one can 
insert each of the following asymptotic expansions 

4> = ^-iC 1 + to + o(C), 

<f> = fa + M + 0(( 2 ), (3.38) 

<f> = M + HC 2 + 0(C 3 ) 

into the Riccati-type equation (|3.16p and, upon comparing coefficients of 
powers of £, which determines the expansion coefficients (frk in (j3.38|) . one 
concludes ()3.33p and ()3.34p . Expansions (13.351) and (I3.36P is an immediate 
consequence of (|3.20p . ()3.33p and (|3.34p . Finally, expression (|3.37p follows 
from (I3TT2]) and (ETI3j> . D 



Lemma 3.4. Suppose that a, f3 satisfy ()3.ip and the pth stationary RT 
system (|2.24[) . Moreover, assume hypothesis (|2.3ip and ()3.2p and let n € Z. 
Let T>jx, fi = {/ti, . . . , fip} be the Dirichlet divisor of degree p associated with 
a,/3, and <j) defined according to (|3.10p and (|3.1ip . that is 

jX{n) = {nj{n), -2Vii(fij(n),n)) G K p , j = 1,... ,p. 

Then T>p i t n t \ is nonspecial for all n 6 Z. 
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Proof. The divisor 'Djxfn) 1S nonspecial if and only if {p,\(n),- • • , p, p (n)} 
contains one pair of {fij,£ij(n)}. Hence, T^p,( n ) is nonspecial as long as the 
projection fij of fij are mutually distinct, fij(n) ^ fJ>k( n ) f° r 3 7^ k. If two 
or more projection coincide for some iiq £ Z, for instance, 

Mji ( n o) = • • • = /tij fc (n ) = /-to, fc > 1, 

then there are two cases in the following associated with hq. 

(i) ^0 ^ {E ,E 1 ,- ■ ■ ,E 2p+ i}; we have Vu(no,n ) / and fi jl (no), ■ ■ • , 
fij k (no) all meet in the same sheet. Hence no special divisor can arise in 
this manner. 

(ii) /io ^ {-^0, Ei,- ■ ■ , £^2p+i}; We assume /^o = Eq without loss of generality. 
One concludes ^12(2, no) = O ((z — Eq) 2 ) and Vn{Eo,no) = 0. Hence 

R2 P +2(z, n ) = - (Hz - l)a + V 12 V+ = O ((A - E ) 2 ) . 



This conclusion contradict with the hypothesis (|3.2p that the curve is non- 
singular. As a result, we have k = 1 and j = are pairwise 
distinct. Then we have completed the proof. □ 

Next, we shall provide an explicit representation of (fi,ipi,a and (3 in 
terms Riemann theta function associated with K, p . 

Let us introduce the holomorphic differentials t]g(P) on K p defined by 

z^ ^ dz 

vi = > £ = i,---,p 

y 

and choose an appropriate fixed homology basis {aj, bj} V j~ 2 on K, p in such a 
way that the intersection matrix of cycles satisfies 

aj ob k = Sj jk , dj o a k = 0, bj o b k = 0, j, k = 1, . . . , r - 2. 

Define an invertible matrix C £ GL(p,C) as follows 

C = (Cj,fc)pxp) = / 

J a k 

c(k) = (ci(k),...,c p (k)), Cj(k) = (C _1 ) i)fc , 
and the normalized holomorphic differentials 



(3.39) 



^2cj(£)r]i, / Uj = 5 jtk , / Uj = T jik , j,k = l,...,p. (3.40) 



18 



One can see that the matrix T = (Tij) pxp is symmetric, and it has a positive- 
definite imaginary part. 

Next, choosing a convenient base point Qq G fC p \ {Poo±, Ph}, the vector 
of Riemann constants Eq q is given by (A. 45) [37], and the Abel maps Aq (') 
and Qq (-) are defined by 

A Qq : JC P ->• J(JC P ) = C p /L p , 



P^A Qo (P) = (A Qojl (P),...,A Qo>p (P)) 

>P rP 



/ UJ!,..., Up) (mod L p ), 
JQo JQo J 



and 



P^«q P)= £ V(P)A Qo (P), 



where L p = {z G C p | z = iV + TM, N, Me V}. 
For brevity, define the function z : fC p x cj 

2(P,Q) = 3q -^q (P) + «q (%), 
PelC p , Q = (Qi,...,Q p ) (3.41) 

here 2(-,Q) is independent of the choice of base point Qo- The Riemann 
theta function 9(z) associated with fC p and the homology is defined by 

= ^exp(27ri < n,z > +iri < n,nT >) , z£ C p , 



ne- 



where < B, C_ >= B_- (7* = ^2j = i BjCj denotes the scalar product in C p . 

(3) 

Let Upp be the normal differential of the third kind holomorphic on 
JC p \{Poo + , Pq} with simple poles at Poo-> Ph an d residues -1,1, respectively. 



1 a p /Cp= JC P x ... x /C p 
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In particular, 



(3) 



y-2V 11 (h-\n)dz + J_^ {z _ Xj)dz 



- i^" 1 

m 2 



z - h- 1 2y 2y ^ 

+|E?=iA i ) + o(C))dc, 



(3.42) 



P -»• Poo-, 



(3.43) 



where the constants {Aj}j =1 6 C are uniquely determined by employing the 
normalization 



(3) 

PhPo, 



o, i = i,...,p. 



The explicit formula (|3.42p and (|3.43|) then indicate the following asymptotic 
expansion near Poo- (using the local coordinate £ = l/z), 



cxp 



Qo 



(3) 

PhPo, 



2p+l 



m=0 
P \ 

+ «E A i) + °(0 , P^Poo- (3.44) 



where Co is an integration constant only depending on K p . Moreover, assume 
r\ E C and \r\\ <min{|i?o|~ \ \E\\ , l-E^I -1 ) - ' ' j l-E^p+il } and abbreviate 



Then 



where 



P = (Eo, E\,- ■ ■ , E2 P +i] 

\ -1/2 



'2p+l 

1(1" ^ntf) 



(3.45) 



vtn=0 



fc=0 
2p+l 



c (P) = 1, h{E) = \ E E «» 



m=0 



&CE) 



j0 f ■ J2p+ 1=0 JoH hi2p+l =fc 



2 2fc (jo0 2 (j 2p+ i!) 5 



E J ...Eg+l, kN, etc. 
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Similarly, 



where 



'2p+l \ 1 / 2 +00 

U(l-E mV )\ = E C *®^ (3.46) 

^m=0 / fc=0 



1 2p+l 

C (P) = 1, Cl (E) = --J2 E m, 



m=0 



Ck(E) = Yl 



(2 J0 )!- • • (2j 2p+1 )!^° . . . ggg 
2^(j !) 2 . . . (j 2p+ i!) 2 (2jo -!)••• (2j 2p+ i - 1) 

fc E N, etc. 



Given these preparations, the theta function representations of 4>,ipi,a, 
and (5 then read as follows. 



Theorem 3.5. Suppose that a, ft satisfy (|3.ip and the pth stationary RT 
system (I2.24p . Moreover, assume hypothesis fl2.31f) and (13. 2D . and let P € 

^p\{-foo±j -Pft} an d (n, no) € Z 2 . Then for each n G Z, £>/i( n ) is nonspecial. 
Moreover, 



0(P,n) 
V>i (P,n,n ) 



where 



and 



0(z(P,A+(n))) 
L (n) — t 777^7 ^ exp 



C(n, ra 



%(P,£(n))) V7Q 
0(z(P,fi(n))) 



(3) 



(3.47) 



0(z(P,fi(no))) 



exp ^(n - n ) ^ 



C{n) 



■ X g(2(Poo-,A("))) 

5(2(^00-, A + (™))) 



(3.48) 
(3.49) 



C(n,n ) 



'IT" 1 C(n' 



n > no, 
1, n = ra , 

{U n jZlC(n')-\ n<n . 



(3.50) 
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The Abel map linearizes the auxiliary divisor T^p,( n ) m the sense that 



"Qo( P A(n)) = ^Qo(%no)) - Ap x A p n){n - n ) 
and a, /3 are the form of 

Ci gU(fio-, £("))) g(l(^oo+^ + (»))) 

£c 0U(^oo-,/£ + (n))) 0(l(Px>+, £("))) 
flU(Poo_, /*"(«))) 0(z(Poc +) A(n))) 



a 



(3.51) 



(3.52) 



and 



P 



(ci/cq)- 



2p+l 



7 E^»-5 R " 1 + 



m=0 



+ 



m 



do;,- 



%(Poo-,/t(n)))+o;) 



LJ=0. 



(3.53) 



Here cq,ci £ C are integration constants. 



Proof. By Lemma 13.44 ^( n ) is nonspecial and hence the theta functions 
defined in this lemma are not identical to zero. Obviously, by (I3.12D and 
Riemann-Roch Theorem 1211 



6(z(P,Kn))) 
LP n) r +( ZTs exp 

v{z\P,H («))) 



(3) 



is holomorphic function on compact Riemann surface K, p and therefore it is 
a constant C{n) related to n. Then <j)(P, n) has the form (j3.47p . Next we 
account for the following Taylor expansion near Poo_ (with local coordinate 

C = l/*), 



9(z(P,(l+(n))) 



^(l(Poo-,/i + (n))) 



0(z(P,p,(n))) C^o eU(Poo_,/i(n))) 



1 + E c i^) 



In 



+ o(c 2 )), 



0U(Poo-,A + N)+ai) 
0(2(-Poo-,/*(n))+w) 

as P ^ Poo_ 



w=0 



(3.54) 
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and hence (|3.44[) and (j3.54p indicate 

0(zfPoo_,/i + (n))) , 
C^o e(z(Poo_,/i(n))) 



2p+l 



1 p fl 



IE* 

m=0 

z(Poo_,/i + (n)) + w)- 



2 



+ 0(0 , 



j'=i 



0(l(Poo_,^(n))) + w) / L=o 

as P ^ Poo_. (3.55) 

A comparison of the coefficients of the asymptotic relations (|3.33|) and (|3.55j) 
then yields the following expressions for C(n) and /3 — ha, 

.! 0U(Px>-,A(n))) 



C(n) 
P — ha 



(3.56) 



m=0 



9 lii / 6>( 1 (P 00 ^,^+(n))+^) 
<9wj V 0(z(P oo _,^(n))) + w) / oj=o. 



(3.57) 



Similarly, one finds the following Taylor expansions near Poo+ (with local 
coordinate £ = l/z), 



0(z(P,£+(n))) 



0(z(P oo+ ,/i+(n))) 



0(z(P,/£(n))) C->o e(z(Poo + ,A(n))) 

9 



<9w. 



In 



i-E c » 

0d(Pco+,A + W) + ^) 



0(l(Poo + ,/i(n)) +w) 



w=0 



+ o(c 2 )), 



as P-)-P 



00+ 



(3.58) 



and 



MP ^ ^fe( P °°+'/f ("))) / 

0(P,n) = ciC n 11 + C2C + O C 

C^o 6»(z(P 00+ ,//(n))) V 

as P -+ Poo+ (3.59) 

where Co, ci are constants arising from the limiting procedure. A comparison 
of (13331 and (Eggp then yields 



ha 



c\C{n 
ci/cq 



e{z{P^+,gn))) 
0(l(Poo-,A(n))) #U(Poo+,A + ("))) 



flU(^oo-,/i + (n))) 0(2(Poo+,A(n))) 



(3.60) 
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which proves (|3,52p . By (|3.57p and (|3.60p . one concludes (|3.53j) . Finally, 
(|3.5ip is the immediate consequence of (|3.12p and Abel's theorem [211 1221 

m\M- □ 

We conclude this section with the trivial case p = excluded thus far. 

Example 3.6. Assume p = 0, P = {z,y) E /Co\{Poo±, Ph} and (n, n ) G 71?. 
Then, 



F + (z,n) 



' -z/2 -ha -(3 + Sx/2 1 \ 

(Hz - l)a+ z/2 - ha + + 8 x /2 J 

/C : P (z,y) = y 2 -(z-E )(z-E 1 ) = 0, E^E^C. 
and a, /3 satisfy 

ha + + ha + fi = -(Eq + E x ) j2, 
(ha+ + (E + £i)/4) (-&*+ + <5i/2) + a+ = -£ £i/4, 

that is, 



P + Pi , 1 {Eq + EA 2 {E -EA 2 E + E l 

a = 2h— ± 2U^h— + 1 - 2~ J ' ^ = ^^ 



Moveover, 



V 12 = 1, K 2 i = (fo- l)a+, 

y n = _ z /2 - ha - p + $i/2, 

V 22 = z/2-Ha + + Ji/2, 

0(P,no) = y/2-{-z/2-hot-p + h/2), 

^(P.n.no) = (y/2-(-z/2-h<x-P + 8 1 /2)} 



n—no 



4 Time-dependent Algebro-geometric Solutions 

In this section we extend the algebro-geometric analysis of Section [3] to the 
time-dependent Ruijsenaars-Toda hierarchy. 

For most of this section we assume the following hypothesis. 
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Hypothesis 4.1. (i) Suppose that a,/3 satisfy 



a(;t),P(;t) G CT, a(n,i)/0, (n,t)eZxR, 
a(n,-),/3(n,-) 6^(1). (4.1) 

(ii) Assume that the hyperelliptic curve IC P satisfies (|2.3ip and f|3 . 2 [) . 

The basic problem in the analysis of algebro-geometric solutions of the 
RT hierarchy consists in solving the time-dependent rth Ruijsenaars-Toda 
flow with initial data a stationary solution of the pth. equation in the hi- 
erarchy. More precisely, given p G No, consider a solution a 1 (n),/3 1 (n) of 
the pth. stationary Ruijsenaars-Toda system s-RT p (a 1 , /3 1 ) = associated 
with ICp and a given set of integration constants {<^}f =0 C C. Next, let 
r G No; we intend to construct solutions a, f3 of the rth Ruijsenaars-Toda 
flow RT r (a,/3) = with a(n, io,r) = a 1 (n), /3(n, io,r) = /3 1 (^) for some 
to r G R. To emphasize that the integration constants in the definitions of 
the stationary and the time-dependent RT equations are independent of each 
other, we indicate this by adding a tilde on all the time-dependent quanti- 
ties. Hence we shall employ the notation V r , RT r , V\i, V12, V21, V22, 9e, fa, $e, 
in order to distinguish them from V p , BT p , Vu, V12, V21, V22, 91, fa, $e, in the 
following. In addition, we will follow a more elaborate notation inspired by 
Hirota's r-function approach and indicate the individual rth RLV flow by a 
separate time variable t r G R. 

The algebro-geometric initial value problem discussed above can be summed 
up in the form of zero-curvature equation 

U tr (z, tr) + U(z, t r )V r (z, t r ) - V+(z, t r )U(z, t r ) = 0, (4.2) 



U(z,t r )V p (z,t r ) - V p + {z,t r )U{z,t r ) = 0. (4.3) 

For further reference, we recall the relevant quantities here (cf. (12 .3p . (I2,10p - 



U= \ „ ^ I, (4-4) 




V{{ V 12 \ „ / V£ V{ 2 



V P =[ " " , V r =\ „ „ , (4.5) 
V y 2i v 22 J V V 22 ) 
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and 



p+i p+i p 

V n = J2g p+1 - jZ j + f p+2 , V12 = Y.h zP+1 ~ j = H( Z -^), 

3=0 3=0 j=l 

fp+1 

V 21 = (hz - l)a + F+ = (hz - l)a+ ]T 

y=o 
p+i 

V22 = -^ffp+l-jZ^, 
3=0 

r+1 r+1 

Fu = £> +1 _ i ^+/ r . +2 , Vi2 = £ (4.6) 

3=0 3=0 

(r+1 
3=0 

r+1 

V22 = -'Y^dr+l-jZ 3 
3=0 



for fixed p e N \{0},r G N . Here {h}Zl,{h}Zl and {Mfij.OtfS 
are denned by (|2.4p - (|2.6|) corresponding to different constants Si and 5e, 
respectively. Explicitly, (|2J| (cf. (f2TT9l) . (IZ20]l ) and (IZ26D are equivalent to 

V n + V u + (z + p)V 12 = 0, (4.7) 
{Hz - l)aV{- 2 -(z + P)V£ - (hz - l)a+V£ + (z + P)V U = 0, (4.8) 
Vn + (z + P) V 12 -V 22 =Q, (4.9) 
«tr + aV£ + (z + P)aV 12 - aV 22 = 0, (4.10) 
Pi, + (hz - l)aV 12 + (z + p)V 22 - (hz - l)a + V+ - (z + p)V u = 0, 

(4-11) 

respectively. In particular, f|2.29j) holds in the present independent setting, 
that is, 

R2 P +2(z) = -V u (z,n,t r ) - Vi 2 (z,n,t r )V 2 i(z,n,t r ), 

= —V n (z,n,t r ) - (hz - l)a + Vi 2 (z,n,t r )V^(z,n,t r ). 

(4-12) 

Here we emphasize that R 2p + 2 is ^-independence (cf. Lemma l4.3j) . 
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As in the stationary context, we introduce 



flj(n,t r ) = (fj,j(n,t r ), -2Vu(fJ>j(n,t r ), n,t r )), j = l,...,p n £ Z, (4.13) 
fXf(n,t r ) = (fj^(n,t r ),2Vii(fif(n,t r ),n,t r )), j = l,...,p neZ, (4.14) 

and note that the regularity assumptions (|4.1|) on a, (3 imply continuity of 
fij with respect to t r £ M. 

In analogy to (|3.10p , (|3. 1 1 [) , one defines the following meromorphic func- 
tion </>(•, n,i r ) on IC P , 



_ y/2 - Vn(z,n,t r ) 
V 12 (z,n,t r ) 

(Hz — l)a + (n, t r )Vi2 (z, n, t r ) 
' y/2 + V n (z,n,t r ) ' 

P = (z,y) G IC P , (n,t r ) eZx 



with divisor (cf>(-, n, t r )) of <f>(-,n, t r ) given by 

(<f>(;n,tr)) = V p ^+(n,t r ) - T> P o-& k n,t r y 



(4.15) 
(4.16) 



(4.17) 



The time-dependent Baker- Akhiezer vector is then defined in terms of <j) by 

' ipi(P,n,n ,t r ,to !r ) \ 
^ 2 (-P,n,no,*r,*0,r) / ' 



(4.18) 



ipi(P,n,n ,t r ,t 0r ) = exp / [Vn(z,n ,s) + V 12 (z,n , s))<f>(P,n , s) J ds 



'Un~=n mn',t r 



n' > no, 



x < 1, 



n = n , 



(4.19) 



Ul.r 



ip 2 (P,n,n ) = exp / ( V\\{z, n , s) + V 12 {z, n , s)(f>(P, n , s) J ds 



x 4>(P, n ,t r ) x < 



T-rn-l ( a(n' ,t r )(hz~l) or I f \ 

lln'=no+l ^ 4>-{P,n',t r ) +z + P\ n i l r) 

1, 

' a(ri ,t r )(hz-l) 



n > no, 
n' = n , 



lln'=n+l ^ <t>-(P,n>,t r ) z~t p\n ,l r ) 
P=(z,y)elC p \{P 00± ,Po}, (n,t r )eZxR. 



n < no, 



(4.20) 
(4.21) 
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One observes that 

ipi(P, n, n ,t r , t , r ) = tpi(P, n , no, U, *o,r)^i n, n , i r , *r) 

^mi 5 »±,Po},(n,no,t r ,io ) r)e2 2 xR 2 . (4.22) 

The following lemma records basic properties of <f> and ^ in analogy to 
the stationary case discussed in Lemma |3. II 

Lemma 4.2. Assume Hypothesis I4.ll and suppose that (|4.2p . (|4.3p hold. 

In addition, let P = (z,y) G /C p \{Poo±, P^}, (n, no, i r , io,r) G x M 2 . Then 
satisfies 

0(P)<T (P) - (z + /3)<T (P) - (fcs - l)a = 0, (4.23) 
tr (P) = (-Kii + F 22 )^ - ^120 2 + (hz - l)a+V+, (4.24) 

^(mn- ^p , (4.25) 

0(P) + ^(P*) = -2^-|j, (4.26) 
W-^) = ^. (4.27) 

Moreover, assuming P = (z, y)\{P 00 ±, P^}, then ^ satisfies 

i/j 2 (P,n,n ,t r ,t ,r) = tpi(P,n,n ,t r ,t 0tr )(j)(P,n,t r ), (4.28) 
17(«)*-(P) = *(P), (4.29) 
V$,(z)*-(P) = (y/2)*-(P), (4.30) 
= K + (z)*(P), (4-31) 
'0i(P,n,n o ,ir,*o,r)'0i(-P* ) n,no,i r ,io,r) = (hz - l) n ~ n °T(a + ,n,n ,t r ) 

x Jl^l) , (4.32) 
Vi 2 {z,n ,t ,r) 

tpi(P, n, n , t r ,t 0ir )ip2(P*,n, n , t r ,t , r ) + ipi(P*,n, n , i r , t ,r)^2(P, n, n , i r , *o,i 

^n-« or , + _ _ ^ Fn(z,n,t r ) 



-2(/iz- l) n - no r(a + ,n,n ,t r 



Vi2(«,n ,to,r) ' 

(4.33) 

-01 (P, n, n , t r ,t 0ir )ip2(P*,n, n , i r , *o,r) - ^i(P*,n, n , t r , t ,r)^2(P, n, n , t r , *o,r) 

= _ l)"-"°r(a+,n,n ,t r ) T . - V - , (4.34) 

Vi 2 (z,n ,to,r) 
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where we used the abbreviation 

f Un~=n f( n '^r), Tl > Tl , 

F(f,n,n ,t r ) = < 1, n = n , (4.35) 

.EKd/K^)- 1 , n<n . 

Moreover, as long as the zeros of fj,j(no,s) of Vi 2 (-,uq, s) are all simple for 
all s € fi, £1 C R, is an open interval, ^ is meromorphic on fc. p \{P QO ±, Pfr} 
for (n,no,t r ) eZx J7 2 . 

Proof. Equations (j4T23"D . (j4T25D - (j4T30|) . (jTOD - fjrai) are proved as in the 
stationary case, see Lemma 13.11 Thus, we turn to the proof of ()4.24p and 
(I4.3ip : Differentiating the Riccati-type equation (|4.23p yields 

<t>t r 4>~ + <N>t T = PtA~ + {z + P)(j)t r + (hz - l)a tr , 

that is, 

- {^ + (^ - (z + P)) S~) <j> tr 

'(hz - l)aV 12 + {z + P)V 22 - (hz - l)a + V+ - (z + p)V 22 ) <fT 
+ (hz - 1) (aV n + (z + P)aV l2 - aV 22 



by using (|4.10p and (|4.1ip . This allows one to calculate the righthand side 

of 



(0- + (^-(z + p)) S')) (Vu<f> + V 12 F - V 21 - V 22 ^ 

= V u <j)<j)~ + Vi 2 <p 2 (p- - V 2l ct)~ - V 22 <\>§- + (hz - l)aV£ 

+ (hz - l)aV^(j)~ - V 21 4> + (z + fi)V 2l - V 22 (hz - l)a 

= Vu {(hz - l)a + (z + p)4r) + v 12 </> {(hz - 1)q + (z + p)<jr) 

- V 2 i<F - V 22 {(hz - l)a + (z + $)()>-) + (hz - l)aV{[ + (hz - VjaV^<jr 

- V 2X <\> + (z + (3)V 2l - V 22 (hz - \)a. 

(4.36) 

Hence, 

(<f>- + (<t>-(z + p)S~)) (-<f> tr + V n <f> + Vi 2 2 - V 21 - V 22 ^) 

= -(z + p)v 22 <\>- + Vu(z + p)tfr + (z + pfv^- ( 4 - 37 ) 

= 0, 
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using (|4.9p . Solving the first-order difference equation (|4.37p then yields 



K + Vn4> + v 12 <t> 2 - v 2 i - v 22 4> 

'lln, =1 B(P,n',t r ), nGN, 
E(P,t r ) x I 1, n = 0, 



(4.38) 



where 



B(P,n',t r 



{P,n',t r ) - (z + /3(n',t r )) 
0-(P,n',t r ) 



(n',t r ) eZx 



and E(-,t r ) is some re- independent meromorphic function on /C p . The asymp- 
totic behavior of (p(P,n,t r ) in (|3.33|) then yields (for t r £ M fixed) 

B(P) = ^C _1 + 0(1) asP^P 00+ . 
(->o to 

Comparing the order of both sides in (|4.38p and taking n > sufficiently 
large, one finds contradiction unless E = 0. This proves (|4.24j) . To prove 
P~3Tj) we rewrite (H~24"l) as 



-Vn + + (z + /3+) y+ - ^120 + 



(fe - 1) a + V& 



-V u + v£ + (z + /3+) - Fi 2 </> + F+<A + - {z + /3 + ) 

^11 + ^12^ + - ^11 - Bl2 ^ 



(4.39) 



Abbreviating 

A(n ,i r 



f().r 



Vii(z,re ,s) + ^2(2;, n o ,s)0(P, n , s) ) ds 



one computes for n > no + 1, 

^i )tr = f exp(A) ' ; 0(n') 



n-l 



n'=no 
n-l 



A t > 1 + exp(A) J{ 4>{n") 
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Vn(z,n ,t r ) + Vi 2 {z,n ,t r )(f)(P,n ,t r )j ipi 
n-l 

+ exp(A) £ (^(^^(nO + V^Cn') 

n'=no 

-Vi 2 (nO0(n')-Fii(n'))0(n') J] <K n ") 

= (^ll + W)Vl- (4-40) 
The case n < no is handled analogously. By (|4.28p and (14.401) . 

= [(-Vu + Vi2)<P - V 12 (f> 2 + (hz - l)a+V&]ih + (f>[Vn^ + V 12 ^ 2 ] 
=V 2 i^i + V 22 ip 2 (4.41) 

Combining (|4T40]) with (JOT]) then yields (j4~3il . 

That ipi(P, n, no, to, t r ) is meromorphic on fC p \{P 00 ±} if V± 2 (-, no, t r ) has 
only simple zeros distinct from P/j is a consequence of (|4.15j) - (j4.17p . (|4.19p . 
03211 and of 

V12 = <9 S In (Vi 2 (2, n ,s)) + 0(1), as z ->• u-,(n , s). □ 
F->-/ij(no,s) 

Next we consider the t r -dependence of Vii, V\ 2 V 2 \- 

Lemma 4.3. Assume Hypothesis 14.11 and suppose that (|4,2p . (j4.3|) hold. 
In addition, let (z,n,t r ) £ C x Z x 1. Then 

Vu ttr =(V n - V 22 )V l2 - 2VuV l2 , (4.42) 

Vl Mr =V 12 V 21 " V 2l V 12 , (4.43) 

V 21 , tr =2V xl V 2l + (V 22 - Vu)V 2l , (4.44) 
In particular, (|4.42j) - (j4.44|) are equivalent to 

V p , tr = [V r ,V p ] (4.45) 
and the spectral curve IC P defined by (I2.3ip and (|4.12p is ^-independent. 
Proof. To prove ()4.42p one first differentiates equation (|4.2T[) 

<^(P)-^P*) = -g^. 

v 12\ z ) 
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The time derivative of eft given in (|4.24p and (|4.26p yields 

^tr(P) ~ A-(P*) =(-Vll + V 2 2)(<P(P) ~ 4>(P1) 

-v 12 (ct>(p)+<p(p*)mp)-HP*)) 

=(-Vn + V 22 )y/V 12 + 2V l2 V liy /V? 2 , 

and hence 

Vi2,t r = (Vu - V 22 )V l2 - 2V n V l2 . 
Similarly, starting from (|4.26j) 

<Pt r (P) + K{P*) = -2%^ + 2 Vl 



V12 V? 2 
yields (|4.43|) . Moreover, 

V 21 , tr =(hz - l)a+y+ + (Hz - l)a+V+ !tr 

= -{hz- l)V+(a + Vu + {z + p + )a + V+ 

- a + V+) + (Hz - l)a + (V+V+ - V+V+ 



-2V+V+) 

= - V 21 V U -{z + /3 + )^2i + V+ - 2V+V 21 
=2V n V 21 + {V 22 - V n )V 21 , 

using gj]) and <^A0\) . Finally, by (jl^) - (Oil) . differentiating (j4TT2l) 

with respect to t r then yields 

R2 P +2,t r = ~ 2V 11 V 11 , tr - Vi2, tr V 2 i - V 12 V 2 i, tr = 0. (4.46) 

□ 

Next we turn to the Dubrovin equation for the time variation of the zeros 
fij of V12 governed by the RT r flow. 

Lemma 4.4. Assume Hypothesis 14.11 and suppose that (|4.2p . (j4.3|) hold on 
Z x with C 1 an open interval. In addition, assume that the zeros 
fij, j = 1, . . . ,p, of Vi 2 remain distinct on ZxZ M . Then {p>j}j=i,...,p defined 
in (I4.13P and ()4.14p satisfy the following first-order system of differential 
equation on Z x 1^, 

p 

H,t r = -Vn{nj)y(fij) JJ (jij - n k )~ l , j = l,...,p, (4.47) 

i=i 
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with 

Ai(n,-) GC°°(X M ,/C P ), j = l,...,p,neZ. 

Proof. It suffices to consider (|4.47j) for fXjj r ■ Using the product representa- 
tions for V\2 in (|4.6|) and employing (|4.13|) and (|4.42|) . one computes 



Mi 

= - 2V n (/ij)Vi 2 (/ij) = y(p,j)Vn{iij), j = l,...,p, 

proving (14.470 . □ 

Since the stationary trace formulas for fi in terms of symmetric functions 
of the zeros /ij of V12 in Lemma [3.21 extend line by line to the corresponding 
time-dependent setting, we next record their independent analogs without 
proof. For simplicity we again confine ourselves to the simplest cases only. 

Lemma 4.5. Assume Hypothesis 14.11 and suppose that (|4.2p . (|4.3p hold. 
Then 

p 

-h(a + a + )-/3 + 5 1 = -Y^ Mj- ( 4 - 4 8) 

i=i 

Next, we turn to the asymptotic expansions of eft and ip\ in a neighbor- 
hood of Poo± an d Pfr. 0.4cm 



Lemma 4.6. Assume Hypothesis 14.11 and suppose that (|4,2p . (j4.3|) hold. 

has the asymptotic behavior 



Moreover, let P = (z,y) £ /C p \{Poo±, P h }, (n,n ,t r ,t r) € Z 2 x M 2 . Then 



MP) iC' + W-h^ + OiC), P^Poo-, , 

C^o \ha + + (h 2 a + a + a-haP)C + 0(C 2 ), P -> P^, 

(4.49) 

,=„ ^T^ + °(C 2 ), P^Pft, C = ^"1A- (4.50) 
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The component ipi of the Baker- Akhiezer vector has asymptotic behavior 
^i(P,n,n A-,io,r) = exp ( ±ht r - t 0>r ) V 5 r+1 . s C is+1) {I + 0(C)) J 



x < 



C°- n (i + o(c)), 



as P -> Pn. 



r(to+) (i + o(0) 

x exp I / J] 5 r+ i_j/ J - +2 (n ,s) ds , 

\j<o,r \j=o y y 

as P -> Poo+, 

(4.51) 

^ 1 (P7i,n ,^,to,r) c = o r(-^ T )C n - no (l + 0(C))ex P n ^ u (ft, n , s)dsj , 

P Pi- 
(4.52) 

Proof. Since by the definition of (ft in (I4.15P and (I4.16P the time parameter t r 
can be viewed as an additional but fixed parameter, the asymptotic behavior 
of (j) remains the same as in Lemma 13.31 Similarly, also the asymptotic 
behavior of ipi(P,n,no,t r ,t r ) is derived in an identical fashion to that in 
Lemma 13.31 This proves (14.51 H for io,r = t r , that is, 

(1 + 0(0), P^Poc-, 



ibi(P,n,na,t r ,t r ) = ( "° , v C = l/z, 

{ ; C->olr (ha+) (1 + 0(0), P^Poo+, 



It remain to investigate 
i/ji(P,n,n ,t r ,t 0)r ) = exp 



(4.53) 



Vn(z,n ,s) + Vi2(z,n o ,s))0(P,n o ,s) ) ds 



(4.54) 

Next, it is convenient to introduce the homogenous representations of Vii, V12 
defined by vanishing all of their integration constants 5i, that is, 
p+1 

Vu = ^ g p+ i-£z e + f r+2 = Vn\s =i,Sj=oj=i,...,p+i, p e N, 

j=0 

P+i 

V12 = U+i-zz 1 = ^i2|5o=i,5 j =oj=i,...,p+i, PGN. (4.55) 

j=0 
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(p) (p) 

In order to avoid confusion about notation, we relabel Vn, Vu by , 

to represent the polynomials associated with the pth stationary Ruijsenaars- 

Toda equation (|2.24p . Then we have 

P+i p+1 
3=0 j=0 

j=0 j=0 



(4.56) 



Focusing on the homogeneous first, one computes as P — > Poo±, 

ytt) , yU) 6 _yU) , f>(i) ~ ^ff 

*n t y 12 <p — v n ^ v i2 _ r (p) 



K 12 



o-) ■ r r u) ^-Vn/v 
vg>/v 



sV & ) + V * MP) 



: E &+i V + /i+a + ( E W ^ 



£=0 



± i r u+i) + //i+a + °(C), ^ 



C = l/*. 



OG- 



2 [0 + 0(0, P^P 

One infers from (|4.56p 

(4.57) 

Insertion of (j4.57p into (|4.54p then proves (|4.5ip . Finally, (|4.52|) follows from 
(tUSD and (530}. □ 

Next, we turn to the principal result of this section, the representation 
of (f>,ipi,a and /3 in terms of the Riemann theta function associated with 
JC p , assuming p € N for the remainder of this section. 

(2) 

Let uj p „ be the normalized differentials of the second kind with a 
unique pole at Poo±, respectively, and principal parts 

w £±, 9 c = ^~ 2 ~ q + d ^ P ~> P -±' C = V*, 9 e No, (4.58) 
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with vanishing a-periods, 



Moreover, we define 

'r+l 



J =1 

,(2) 



and corresponding vector of 6-periods of /(2ni) is then denoted by 



bj 



Finally, we abbreviate 
= lim 

and 



(/>4t>-<-<- +i) ) 



P-+P00- 



Qo 



Theorem 4.7. Assume Hypothesis 14.11 and suppose that ()4.2|) . (|4.3|) hold. 

Moreover, let P = (z,y) € /Cp\{Poo ± , P R }, (n, n , t r , to,r) G x ]R 2 , Then 
for each (n, i r ) £ZxR, T^jx{n,t r ) ls nonspecial. Moreover, 



4>(P,n) =C{n,t r ) 



6(z(P,fi + (n,t r ))) 



%(P^(n,t r ))) eX </ Qo ^-)' (4 ' 60) 
ipi{P,n,n ,t r ,to tr ) =C{n,no,t r ,t , r ) —. , n —, 7 — TTT ex P ( n ~ n o) / 

-(*r-V) / P ^ 2) )> ( 4 - 61 ) 



PftPx>- 



where 

~, , v -1 9(gPoo-4^htr))) 



3G 



and 

C(n, n , t r ,t 0>r ) = exp ((t r - t ,r)f^°~ ~(n- n )uj ^ 

6»(z(Poo_,/i(no,to,r))) 



(4.63) 



6(z(Poo-,[l(n,t r ))) 
The Abel map linearizes the auxiliary divisor "Dfi( n ,t r ) in the sense that 

~ ( 2 ) 

OQ (%n,t r )) = a Qo {V^ n0ttor) )-A Pao _(P h )(n-n )-U r (t r -t , r ) (4.64) 
and a, [5 are the form of 

+ , ci g(l(Poo-,Mn,t r ))) g(l(Poo+^ + (n,t r ))) 

" (n ' r) hcoe(z(Poo-,p. + (n,t r ))) 9{z{P^Kn,t r ))) l ' j 



and 



0/ . x / / 0{z{Poo-,V (n,t r ))) 0(^00+, p,(n,t r ))) 
p{n,t r ) = (ci/co) 



0U(Poo- , /2(n, i r ))) 0(2(P oo+) pT (n, tr))) 

1 2p+l p 
m=0 j'=l 



9 , /^(l(Poo-,/i + (n,tr))+w) 

+ — l n 



3wj V ^(^(Poo-, fi(n,t r ))) + w) 
Here Cq,Ci € C are integration constants. 



(4.66) 

uj=0. 



Proof. As in Theorem !3.5l one concludes that 4>(P, n, t r ) has the form (j4.60|) 
and that for to, r = t r , ipi(P, n, no,t r , to,r) is of the form 

ipi(P,n,no,t r ,t r ) = C(n,no,tr,t r )- 



6(z(P,tl(n Q ,t r ))) 
p 



exp ({n - n ) J 4^+) • 



'Qo 

To discuss ipi(P,n,no,t r ,to tr ) we recall (|4.22[> . that is, 

V>i(P,n,n , Wo,r) = Vi(- p 5 «o,n ,ir,*o,r)V'i(- p , n,n ,t r ,t r ), (4.67) 
and hence remaining to be studied is 

i>i(P,no,no,t r ,t Qjr ) = exp ( / V n (z,n , s) + V 12 (z,n , s)<j}(P,n , s) J . 

(4.68) 
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Introducing ip\(P) on JC P \{P 00 ±} by 



C(riQ,no,t r ,to jr ] 



x exp 



e(z(P,jj(no,t r ))) 
d(z(P,fi,(n ,t 0j r))) 
■p 



-(tr-t ,r) / ^ 



Qo 



(4.69) 



we intend to prove that 



ipi(P,n ,no,t r ,t ,r) = i>i(P,n ,t r ,t 0i r) 
P G /C p \{Poo ± }, n G Z, t r ,t , r G R : 



(4.70) 



for an appropriate choice of the normalization constant C(no,nQ,t r ,tQ r ) in 
P~69~|) . We start by noting that a comparison of (|43T]l . (Ii32j) . (|459|) . (^6ll 
shows that tp\ and ipi have the same essential singularities at -Pooi- Thus, 
we turn to the local behavior of ipi and By (14, 69ft ifti has zeros and 
poles at fi(no,t r ) and /i(no,io,r)- Similarly, by (|4.68|) . has zeros and poles 
only at poles of cj>(P,no,s), s G [to,r> t r ](resp., s G [t r) io,r])- In the following 
we temporarily restrict to,r and t r to a sufficiently small nonempty interval 
ICR and pick no G Z such that for all s G I,/ij(no,s) ^ Ufc(no,s) for all 
j ^ k,j,k = 1, . . . ,p. One computes 



if>i(P,nQ,n>o,t r ,t ,r) 

exp / ds [V n {z,n ,s) + Vi 2 {z,n ,s) 



to,r 
( 



P— tp,j(no,s) 



exp 



ds 



t().r 



Vu(z,n , s) 
Vi2 («j(n , s)) y{fij{n , s)) 



+ 0(1) 



(z - fij(n ,s)) II («i(«o,s) - u fc (n ,s)) 
fe=i 



P^>-jlj(n ,s) 



P-¥p,j(no,s) 



exp 



exp 



to 



to 



ds 



ds 



z - Hj(n ,s) 
d 



+ 0(1) 



ds 



In (nj(n ,s) - z) + O(l) 



(4.71) 



Restricting P to a sufficiently small neighborhood Uj(no) of {/ij(no,s) G 
JC p \s G [io,r,*r] -0 such that (J,k(no,s) £ Uj(no) for all s G [io,r,*r] ^ ^ and 
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all k G {1, . . . ,p}\{j}, and (13771)1 imply 

^i(P,n ,n ,t r ,t 0j r) 

(fij(n ,t r ) - z)0(l) , as P -> jj,j(n ,t r ) / Aj( n o, *0,r), 

0(1), as P fLj(no,t r ) = /xj(n ,io,r), 

(^j(n ,to,r) - (1) , as P jj,j(n ,t 0>r ) ^ p,j(n ,t r ), 

P = (z,y)elC p , (4.72) 

with O(l) ^ 0. Thus and have the same local behavior and iden- 
tical essential singularities at Px,±. Hence an d coincide up to a 
multiple constant (which may depend on no,t r ,to tr ). By continuity with 
respect to divisors this extend to all no £ Z since by hypothesis T>^ n ^ 
remain nonspecial for all (n,s) € Z x R. Moreover, since by (|4.68p . for 
fixed P and no, ^>i(P, no, no, io,r) is entire in t r (and this argument is sym- 
metric in t r and io,r), (|4.70p holds for all i r ,to,r- € M(for an appropriate 
choice of C(no, no, i r j to,r))- Together with (I4.67p . this proves (|4.6ip for all 
(n,t r ), (n ,i ,r) eZxR. 

To determine the constant C(n, no, t r , io,r) one compares the asymptotic 
expansions of Y>i(P, n, no, t r , io,r-) fo r P - ^ Poo- m (|4.5ip and (|4.6ip 



C(n, no, i r , *o,i 



Finally, (|4.64j) follows from 
d 



exp Ut r — io,r)^£° — (n — no)^o 

^(^(Poo-, A( re 0,*0,r))) 

^ — 

9(z(Poc-,p.{n,t r ))) 



dt, 



^Qo^( P A(n,tr)) 
p,j(n,t r ) 

Qo 



J = l 



0J£ 



p / p 

j=i \fc=i 

V ( V 

i=i \fc=i 



.fc-i 



y(fij(n,t r )) 



V 



Fi 2 (Mj(n,t r ))y(Ai(ra,t r )) 0^ 

Vl2(Mj(ra,tr)) 



"fcj 



fc=l 
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k-l 



^ ^ Q(A;)5 r _ s c fc+s _p( ; E) 



fc=l i 

p r 



r+1 
s=0 



1-s 



. t=max{0,s— p} 



fe=l s=0 
(2) 



where we use the interpolation representation of V12 in appendix A (cf. 
g29])) and 



c 



p-j 



ifiS(i 

00 p 



g=0 fc=l 



(2) 



1 



r ^ 2vri 



6j 



1 rl r+1 

S=l 



(2) 

Poo+,S-l 



(2) 

Poo-,8-1 



r+1 



: ^ £ r +l-8 ^ Ci(k)c k _ p+S (E). 



□ 



8=1 



fc=l 



Appendix A: The Lagrange Interpolation Represen- 
tation of Vi2(/ij(n, t r )) 

Introducing the notation in [371 [38], 

*fc(M) = (-1)* •••,/%, 

£G5 fc 

= {I = (4, • • • ,4) e N fc |4 < . . . < 4 < P}, k = l,...p, 
= (-l) fc X •••,/%, 

r W = {£ = (£ 1: . . . , 4) G N fc |4 < . . . < 4 < p 4, ^i}, 

fc = i,...,p-i, i = i,...,p. 
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and the formula 

k 



^* fc _/(M)^ = $^(M), k = Q,...,n, j = l,...,n, (4.73) 



one finds 



p+i p p 

Vi2(2) = E /p+1-52* = II (* - Hj) = e vp-tin)** 

s=0 j=l 1=0 

and 

ft = *t-i(ii), e = 0,---,p + l. (define *_i(m) = 0) 
In the case r < p, 



r+l 

r+l / min{r+l— s,p+l} 



^12=X)/r+l- 
r+l / min 

= E E c r+1 _ s _ fc GE)/ fe ]z s 

s=0 \ fc=0 / 

r+l /r+l-s \ 

= E E Zr+l—k(M>fk K 
s=0 V fc=0 / 
r+l /r+l-s \ 

= E E Cr+l-— Jfc(£)**-l(£) 2 S 
s=0 V fc=0 / 
r+l r+l— s 

= E^^) E ^r-s-tiriz* 
s=0 t=0 
r r—s 

= 5>(£)X>r-,_ t 0i)s*. 



s=0 t=0 
Using (|4.73p . we have 



*Wj) = Yl E *r- S -t(/f)M 

s=0 t=0 

s=0 



(4.74) 



(4.75) 
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In the case r > p, 

V 12 (z) 
r+l 

= E fr+l~sZ S 
s=0 

r+l / min{r+l-s,p+l} \ 

= E E c r+1 - s - k (E)f k z s 

s=0 \ k=0 J 

r—pp+1 r+l r+l— s 

= ^2^2c r+1 - s - k {E)^ k ^ 1 (fl)Z S + E E Cr+l^fc(^)*fe-l(^)2 S 
s=0 fc=0 s=r-p+l k=0 

r—pp+1 r+l p+1 

= ^2^2c r+1 . s . k {E)^ k _ 1 ( [ i)z s + E 5^c r+1 _,_ fc CE)* fc _i(/i)js a 

s=0 fc=0 s=r-p+l fc=0 

p+1 r+l 

= E E Cr+l-s-fc(^)^fc-l(/f)^ 
fc=0 s=0 
r+l p+1 

= E E <w-s-fc(£)*fc-i<>)^ 

s=0 fc=0 
r+l p+1 

=EE^cs)**-i« r+1 — * 

s=0 fc=0 

r— p /p+1 \ 

s=0 \fc=0 / 

r+l /p+1 \ 

+ E p.® E**-^ 2 ^ 1- *"* 

s=r-p+l \fc=0 / 

r— p r+l /r+l— s \ 

=E^)cm*)K~ p_s + E E ^-^> r+1 - s - k )- 

s=0 s=r-p+l \ fc=0 / 

(4.76) 

Then one finds 

r+l /r+l— s \ 

^2(^0= e e **-i^r l- * -fc 

s=r-p+l V fc=0 / 

r+l 

= E m)^l(a)- (4-77) 

s=r— p+1 
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Combining (|4"775j) with fljTTD yields 

Vi2{z)= E umiilsM = E ^)*2>) (4.7s) 

s=max{0,r— s=max{0,r— p} 

and hence 

r+l 
s=0 

r+l / S \ ( 4 ' 79 ) 

s=0 \t=max{0,s-p} / 

Appendix B: Asymptotic Spectral Parameter Ex- 
pansions 

Next, we turn to asymptotic expansions of various quantities in the case 
of the Ruijsenaars-Toda Hierarchy. Consider a fundamental system of so- 
lutions *±0, •) = (ipi,±(z, -),tp2,±( z ^ ')) T of U(z)^±(z) = * ± 0) for z E C 
(or in some subdomain of C), with U given by (|2.3p . such that 

det(tt_(*), *+(*)) ^0. 

Introducing 

0± = feM zec> nGN) (4 . 80) 

then cp± satisfy the Riccati-type equation 

<f>±fe-(z + P)<f>i-(hz-l)a = 1 (4.81) 
and one introduces in addition, 

f= , 1 t • (4.83) 

<P+ - <P- 

Using the Riccati-type equation (|4.8ip and its consequences, 
(f> + (f>+ - -{z + p){cj>- - <j)Z) = 0, 

4>Z(f>+4>+ - 4>-4>Z(f>+ = a{hz - i)(cf)Z - <j>+), 
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one derives the identities 



5 + g- + (z + /3)f = 0, (4.84) 

(Hz - l)af~ -(z + /3)g~ - (Hz - l)a+f+ + (z + /3)g = 0, (4.85) 

g 2 -a + {hz-l)ff + z = l/A, (4.86) 

f_ 1 _ ^ = ( 487 ) 
<t>+~4>- a(hz-l)((f)Z - (/>+)' 

Moreover, (14,84j) -( j4.87p also permit one to derive nonlinear difference equa- 
tions for f and g separately, and one obtains 

(z + f3) 2 - (g + g+)( + Q -)(hz - l)a+ (4.88) 

= \(z + P)\ 

({z + /3) 2 f + (Hz - l)af- - (Hz - l)a+f+) 2 

- 4ff + a + (hz -l)(z + f3) 2 = {z + (3) 2 . (4.89) 

Theorem B.l Assume (|2.24[> . s-RT p (a,/3) = 0, and suppose P = (z,y) G 
fC p \{P 00 ±}. Then f, g has the following convergent expansions as \z\ — > oo, 



±X>C', C = V*, 



f^E/^, C = iA, 

1 ^ 1 £=0 



(4.90) 



and simultaneously as P — > Poo±, 

oo 

V u /y = ±J29iC e , C^0, C = V2, 

oo 



(4.91) 



Moreover, one infers for the E'm-dependent summation constants 6e,£ = 
0, . . . ,p, in Vij(i,j = 1,2) that 

5 t = ct(E) £ = 0,...,p (4.92) 
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and 



l 

h=Y, c t-k(R)L £ = 0,...,p + l, (4.93) 

fc=0 

min{^,p+l} 

h= E ct-k{E)fk 

k=0 
min{f ,p} 

= E & £-k(E)fk, (4.94) 

k=0 

Proof. Identifying 

*+(*,-) with #(P,-,0) and *_(*,-) with *(P*,-,0), (4.95) 

and similarly, identifying 

0+(z,-) with 0(P,-) and <j)-(z,-) with </>(P*,-), (4.96) 

a comparison of (|4.80p - (|4.89p and the result of Lemma 13.11 and 13.31 shows 
that we may also identify 

Vn V12 
q with ± , f with ± , 

y y 

The sign depending on whether P tends to Poo±- Hence we are only to 
investigate the asymptotic expansions of V\\jy and V^/y- Dividing Vn and 
V12 by y, one obtains 

VU{Z) - (f:*(®z-*) (^Stz-A = f>-', (4.97) 



Vi 2 (z 



^fc=o / v=o / e=o 

' oo \ /p+1 \ oo 

X>csK fe E^~' = E^~ Z > ( 4 - 98 ) 



for some coefficients fa to be determined next. Dividing (|2.19p and (I2.20P 
by y and inserting the expansions (|4.97p and (|4.98p into the resulting equa- 
tion then yield the recursion relations (|2.4p - (|2.6p (with fa replaced by fa). 
Moreover, plugging (g^Tj) and flOSD into (f2732|) and (g^D then yields 50 = 
-5 = 5o, 5i = ha + = 9i, fo = = /o, /1 = 1 = /1. Then we can induc- 
tively to show that g e = M e (go, ■ ■ ■ ,9i-i) = k, fa = ^e{fo, fe-i) = fa, 
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where 5Dt£,9fy are the polynomials in (50, ■ ■ ■ , §£-1) and (/o, . . . , fi-i), re- 
spectively. This implies (|4.90j) and (|4.9ip . A comparison of coefficients in 
(|4.98|) then proves (|4.94|) . Next, multiplying f|3.45j) and (|3.46[) . a comparison 
of coefficients of rj~ k yields 

Hence one computes 



k + 



I i /min{fc,p+l} 

fc=0 m=0 \ s=0 



k=0 \s=0 J 

I k 

k=0 s=0 

I I 

J2J2 C£ -k(^> d k-s{E)f s 

k i 
s=0 k=0 

s=0 \k=s J 



=ft, £ = 0,...,p+l. 
Hence one obtains ([HMD and tllus (IHS2D- □ 
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